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ABSTRACT 



The Naval Research Laboratory's Far Ultraviolet Cameras 
Experiment is part of an ongoing project to develop a means 
for real time space based monitoring of the upper atmosphere. 
The purpose of this paper is to develop a method by which 
useable data can be extracted from the images obtained by the 
cameras. The required output is an intensity profile of the 
light in the band of the camera as a function of the observer 
zenith angle. From this data future analysis can be done to 
determine the spatial distribution and density of the 
constituent species of the upper atmosphere. The long term 
application of this project is to develop a method for using 
space based imaging of the ionosphere to provide real time 
radio frequency propagation analysis. 
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INTRODUCTION 



A. BACKGROUND 

The purpose of this thesis is to develop a means for 
analyzing images obtained by two far ultraviolet cameras 
developed by the Naval Research Laboratory. The Far 
Ultraviolet Cameras Experiment was part of the Air Force 
Program 675 (AFP-675) payload flown aboard space shuttle 
mission STS-39 between 28 April and 6 May 1991. AFP-675 
consisted of the following instruments: 

1) CIRRIS - Cryogenic Infrared Radiance Instrumentation 
for Shuttle. This instrument measured spatial, spectral, 
and temporal variations in the Earth's limb and aurora 
between 2.5 and 25ixm. 

2) FAR UV - Far Ultraviolet Cameras. Designed to image 
near Earth and celestial ultraviolet sources, including 
airglow, aurora, stars, nebulae, shuttle glow, shuttle 
thruster firings, and shuttle HjO discharges. This paper 
will address the airglow portion of this experiment. 
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3) URA - Uniformly Redundant Array. Demonstrated coded 
aperture imaging of x-ray emissions from celestial 
targets. 

4) HUP - Horizon Ultraviolet Program. Measures spatial 
and spectral ultraviolet characterisitics of the Earth's 
horizon. 

5) QINMS - Quadrupole Ion and Neutral Mass Spectrometer. 
Monitors contamination in the shuttle environment and 
provides input to CIRRIS. (STS-39 PID) 

B. FAR ULTRAVIOLET CAMERAS EXPERIMENT 

The objectives of the Far Ultraviolet Cameras Experiment 
included; far ultraviolet (1050-2000A) measurements of the 
Earth's upper atmosphere, celestial bodies including diffuse 
background sources, and the shuttle environment. The shuttle 
enviroment study included shuttle glow, and thruster firings. 
Upper atmosphere observations include night horizon airglow 
and the auroras. The purpose of these observations was to 
demonstrate a means for remote sensing of atomic oxygen and NO 
density versus altitude in the upper atmosphere. The night 
airglow observations were specificaly designed to demonstrate 
the ability to conduct remote sensing of the night ionosphere. 
The results of this experiment will be used to support future 
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space based real time remote sensing of the night ionosphere. 
This will be attempted by the Global Imaging Monitor of the 
Ionosphere (GIMI) project scheduled to be deployed aboard the 
ARGOS P-91 satellite. GIMI will use similar cameras but will 
utilize CCD detectors instead of electrographic film to 
provide a real time digitized RF output. This output will be 
processed at a ground station to determine an altitude 
profile of the atmosphere. The work done in this thesis will 
provide the basis for deriving altitude and intensity from the 
photographs and will provide the background material necessary 
for further research on this project. The methods used to 
analyze the images from the AFP-675 project should be similar 
to the methods required to analyze the GIMI output; therefore, 
the methods covered herein should support GIMI as well as AFP- 
675. (Carruthers, 1992) 

C. THESIS OUTLINE AMD OBJECTIVES 

In addition to the purpose stated above, this paper is 
intended to provide necessary background information for 
follow on work. The paper is organized as follows. 

Chapter I: This chapter provides general background 
information on the Far Ultraviolet Cameras Experiment and 
future applications of work associated with the project. 

Chapter II: This chapter provides background information 
on the structure of the ionosphere, chemical and dynamical 
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processes in the ionosphere, and a brief introduction to the 
night airglow phonemena. Topics covered include diffusion 
mechanism, density scale height, ionization rates, and 
photochemical equilibrium. 

Chapter III: Provides a detailed description of the Far 
UV Cameras and the methods used to reduce the image data into 
data files. 

Chapter IV: Covers data analysis performed as part of 
this thesis. Includes a description of the means by which 
shuttle position and attitude were determined and details the 
programs used to interpret the image data files. This section 
includes graphs of intensity versus observation zenith angle 
obtained from the images. 

Chapter V: Contains a summary of findings and topics for 
future NPS students to research. 
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II. BACKGROUND 



A. THE UPPER ATMOSPHERE 

1. General Description 

a. Temperature Structure 

The upper atmosphere can be described in many 
ways. The most common means of describing the layers of the 
atmosphere is by the temperature profile. Figure (1) shows a 
theoretical plot of temperature versus altitude and identifies 
the various layers. 




TEMPERATURE K 



Figure 1: Temperature Profile of atmospheric 

layers . 
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In the thermosphere the sun's radiation in the vacuum 
ultraviolet region is absorbed by atomic and molecular 
nitrogen and oxygen. This absorption generates heat, as 
altitude decreases less vacuum ultraviolet radiation reaches 
the atmosphere resulting in lower temperatures. This explains 
the positive temperature gradient in the thermosphere. The 
temperature in the thermosphere reaches a minimum at about 85 
km (The mesopause) . Below this point ozone becomes an 
important constituent and it absorbs longer solar radiation 
wavelengths (2000 -300oA) . From 85 km to about 50 km the 
increasing density of ozone causes a temperature increase with 
decreasing altitude. Also absorbed by ozone is the Lyman-a 
radiation at 1216A which passes through the nitrogen and 
atomic oxygen to an altitude of about 65 km. At 50 km enough 
of the 2000 - 300oA and Lyman-a radiation has been removed to 
allow a positive temperature gradient to form again even 
though the density of ozone continues to increase. 

This thesis will deal primarily with the 
atmosphere above 100 km. Therefore our discussions will be 
limited to the thermosphere. A profile of thermospheric 
temperature is shown in Figure (2) . This profile was 
calculated using the equation: 

T-T.- (r„-Ti2o) (1) 
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where T, is the temperature at the outer edge of the atmosphere 
and T, 2 o is the temperature at altitude 120 km, z is altitude, 
and s is .0255 km'’. The figure shows temperature calculated 
in 2km increments from 120 to 300 km. T„ can vary from lOOOK 
and 1800K depending on the solar cycle. 1500K is used in this 
calculation. A value of 200K is used for T^jq. 




200 400 BOO too 1D00 1200 1400 1B0O 



Temperotire in Kelvin 

Figure 2: Temperature vs Altitude from 120-300 km. 

Too = 1500 K, T^20 = 200 K 

b. Diffusion Structure 

Another means of describing the atmosphere is by 
the mechanism of diffusion which dominates at a given 
altitude. Their are two basic mechanisms of diffusion for the 
neutral atmosphere, molecular diffusion, and eddy diffusion. 
Eddy diffusion is the turbulent mixing of parcels of air in a 
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given region. This is the dominant method in the troposphere, 
stratosphere, and mesosphere where the atmosphere is more 
dense and acts like a fluid. In the troposphere and 
mesosphere strong negative temperature gradients allow the 
easy rise of warm air and sinking of cooler air. 
Alternatively, the dominant diffusion method in the 
thermosphere is molecular diffusion. In this region the 
strong positive temperature gradient inhibits the vertical 
mixing of the atmosphere. The one exception to this is the 
region between 85 km and 120 km where atmospheric tides create 
strong vertical wind shears in the prevailing horizontal 
winds. The stratosphere has a temperature gradient similar to 
that of the thermosphere; here the molecular diffusion rate 
approaches the eddy diffusion rate but eddy diffusion still 
dominates. 

As a result, the atmosphere can be divided into 
two diffusion regions; the homosphere and the heterosphere. 
The homosphere is the region below 90 km and has a nearly 
constant mean molecular mass. The region above 115 km is 
called the heterosphere and is characterized by the mean 
molecular mass decreasing with altitude due to varying scale 
heights. (Scale height will be covered in Chapter III.) 
Figure (3) gives a general indication of diffusion types and 
rates versus altitude, and illustrates the two regions. The 
figure illustrates the crossover point at 115 km and shows the 
variable nature of the eddy diffusion term. Between 90 and 
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(Heicklen, pg 4) 



115 kin molecular diffusion and eddy diffusion are comparable 
as illustrated in Figure (3) . 
c. The Ionosphere 

The ionosphere is so named because of the 
ionization of its constituent elements by solar radiation. 
The ionosphere is described by regions of varying electron 
density. The ionoshere is divided into four distinct regions; 
the D region which is lowest in altitude, the E region, the 
region, and the region which is the highest. The D region 
extends from about 50 km to 90 km blending in with the lower 
region of the E region with no sharp boundary. The D region 
experiences total solar control. It builds rapidly in the 
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morning, maximizes near local noon, and rapidly diminishes at 
sunset. The E region extends from the top of the D region to 
the bottom of the F region, approximately 90 - 130 km. The 
electron density of the E region peaks at about 110 km. The 
E region is very predictable reaching a maximum slightly later 
in the day than the D region. The E region does not disappear 
at night but does reach a minimum density. The F region is 
the highest ionospheric region. The F region extends upward 
from about 130 km and has no clearly defined upper limit. 
During daylight the F region splits into two distinct layers, 
F, and F 2 . The F^ layer is the lower and less dense layer. It 
is controlled directly by solar activity and disappears at 
night blending into a single F layer. The F 2 region is the 
highest and most dense layer with a daytime peak altitude of 
approximately 300 km. The F 2 region is the most variable of 
all the ionospheric regions. The F region reaches maximum 
density a few hours later than the D or E regions and does not 
disappear at night. This is due to the substantially lower 
loss rate of the constituent ions in the F region compared to 
those in the D and E regions. A graph of electron densities 
is given in Figure (4) . 
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Densities . 



2. The Mass Continuity Equation 

There are several mechanisms for the production and 
loss of the constituent species of the ionosphere; diffusion, 
chemical reactions, and photoionization are the major 
mechanisms. The mass continuity equation describes production 
and loss involving these mechanisms. The one-dimensional 
approximation to the mass continuity equation is: 



dn „ r 5T 
dt oz 



( 2 ) 
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where dn/dt is the rate of change in number density, P is the 
chemical production rate, L is the chemical loss rate, and n 
is the density of the specific species, and 8T/dz represents 
loss due to diffusion. Each term will be discussed in detail 
beginning with diffusion. 

a. Diffusive Equilibrium 

The major neutral species (0, 0^, N 2 ) are primarily 
affected by diffusion. Chemical production and loss also 
occurs with these species; however, their rates are so small 
in comparison to the rate of diffusion that they can be 
ignored (P-nL<8T/dz) . A state of diffusive equilibrium exists 
when diffusion dominates and dn/dt = 0. This results in 
equation (2) reducing to; 

0--|t (3) 

dz 

From this equation it can be shown: 



n(z) -ii^exp ^ 

where n(z) is the number density of the species and has units 
of cm'^. Z is the altitude at which the density is calculated 
and is the lower altitude limit (120 km for this paper). 
H* is the density scale height and is given by; 
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1 



( 5 ) 






kT T dz 



where m is the atomic mass of the atom or molecule, g is the 
acceleration due to gravity, k is Boltzmann's constant, T is 
temperature (Kelvin) , and dT/dZ is the temperature gradient. 
Scale height is the difference in altitude that it takes the 
density of a species to drop by 1/e. 



b. Photochemical Equilibrium 

For the major ionic (N 2 *, O^* , O* , NO*) , and the 
minor neutral (NO, N(^S), and N(^D)) species in the ionosphere, 
chemical reaction rates greatly exceed diffusion rates; 
therefore, dT/dz <P-nL and diffusion is neglected. A state of 
photochemical equilibrium exists when chemistry dominates and 
dn/dt=0. Thus chemical production and loss rates are equal 
(P-nL=0) and equation (2) becomes: 




( 6 ) 



c. Photoionization 

Photoionization by solar radiation is the primary 
means of ion production in the ionosphere. The highest rate 
of ion production results from the excitation of atoms and 
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molecules by solar radiation. The ionization rate (Q) of a 
given species is; 



(?(z,0) “n(z) o_j(X) J;^(z,0) AA 



(7) 



where n(z) is the number density given above, is the 

ionization cross section, and 1^^ is the intensity of the solar 
radiation, at the wavelength A., at the altitude z, when the 
solar zenith angle is 0. Each of these items is described in 
detail in the following paragraphs. 

The intensity can be calculated using the 
following equation: 

Ji(z,0) (8) 
where Ij (°°) is the intensity at the top of the atmosphere and 
t(z,0) is the optical depth. is tabulated. 

The optical depth can be calculated using the 
following equation: 

N{z,Q) (9) 

where o^{k) is the species absorption cross section at the 
wavelength X and N(z,0) is the slant column density. The 
optical depth is a measure of the attenuation of solar 
irradiance by the Earth's atmosphere. Optical depth is 
dependent on the specif ic^Xxavelength, the absorption cross 
section of the given species, the altitude dependent 
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concentration of the given species, and the solar zenith 
angle. Where optical depth reaches unity (t= 1), a noticeable 
drop in intensity (Ij^(Z)) occurs. This drop is illustrated in 
Figure (5) . 




OPTICAL DEPTH 

Figure 5: Optical Depth vs Intensity. 



The column density, N(z,0), of the species and can be 
calculated from the number density using the following 
equation; 



N{z, 0) 




dz 

COS0 



( 10 ) 
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where 0 is the angle between the direction of the solar 
radiation and the vertical axis, generally called the solar 
zenith angle. 



3. Chemical Production and Loss in the Ionosphere 

Using the ionization rate described above, the mass 
continuity equation, and using tabulated equations for the 
production and loss of ions in the ionosphere, it is possible 
to calculate the density of each ionic species. The constants 
associated with these equations for production and loss are 
given in Appendix A. Equations for determining the density of 
the major (N 2 *, 02"", 0*, NO*) and the minor ( NO, N(*S), and 
N(^D)) constituents are given below. The reaction 
coefficients are listed in Appendix A. 



The following reactions produce 0*: 



0+h\ - o*+e 



( 11 ) 



N 2 +O - 0"+i\^2 



( 12 ) 



The following equations account for the loss of O*: 



^ N0*+N(^S) 



(13) 
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0* + 02 -* O 2 +O (14) 

Using equation (8) and reactions (11) through (14), the 
density of O* can be calculated using the following equation: 

[0*1 - (15) 

where Q is the ionization rate discussed above. The subscript 
following the reaction coefficient K is the reaction number 
listed in Appendix A, and [X] is the number density of species 
X. Similar calculations can be done for NO*, NO, N(^S), N(^D), 
N 2 *, and O 2 *. 

NO* is produced by the following reactions: 



N2+O -* NO^+N{^D) (16) 

OVW 2 NO*+N(*S) (17) 

O2+NO - NO*+N(^S) (18) 

02+A^(^5) - N0*-*-0 ( 18 ) 



The following reactions lead to the loss of NO*: 

NO* + e - 0+N{*S) (20) 
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iVO* + e - 0+N(^D) 



( 21 ) 



Applying reactions (16) through (21) , the density of can 
be calculated from the following equation: 

r ^ 17 a [^2] [O] +^^20 [O*] [JVj] 10;] [//O] +X23 [q^] [N(*S] ^ 

' tel 

The following reactions account for NO production: 



N{^S) +O 2 ^ NO+0 ( 23 ) 

N{^D) +O 2 ^ NO+O ( 24 ) 

The reactions for NO loss are: 

Oz-^-NO - NO* -^02 ( 25 ) 

Ni^S) +NO ^ N 2 +O ( 26 ) 

N(^D) +N0 ^ N 2 + O ( 27 ) 



Combining reactions (23) through (27) yields the following 
calculation for the number density of NO: 



K26 [N(^S) ] [O2] +K2^ iNj^D) ] [O2] 
K 22 [O;] +J^33 [N{^S) ] +^^34 [N{^D) 



( 28 ) 
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The reactions for N(*S) production are: 



0*+N2 - NO*+N{^S) 


( 29 ) 


NO^e - 0^-N(*S) 


( 30 ) 


N(^D) +e - N(*S) +e 


( 31 ) 


N(^D) +0 - N(^S) +0 


( 32 ) 



N(^S) loss results from the following reactions: 

O - MP" + 0 (33) 



N(^S) +O 2 - WO+0 


( 34 ) 


N(*S) +N0 - N 2 + O 


( 35 ) 



Combining these equations in the form of equation (8) yields; 






[0*] [N^] [ej -*-^30 ] [e] [Nj^D) ] [O] 

[O;] [O^] [NO] 



( 36 ) 



N(^D) is produced by the following reactions: 
Nz+O - NO^+N{^D) 

Nz+e - N(^D) +N(^D) 



( 37 ) 

( 38 ) 
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NO* + e - 0+N{^D) 



( 39 ) 



The following reactions result in N(^D)loss: 



N(^D) +O 2 


- NO+0 


( 40 ) 


Ni^D) +e - 


N{*S) +e 


( 41 ) 


N(^D) +0 - 


N{*S) +0 


( 42 ) 


Ni^D) +NO 


- 


( 43 ) 


Combining equations (37) through (43) and solving for 


[N(2D) ] 



yields: 



[N{^D) ] 



^I 7 a [^2 ] [O] [N^] [e] +7^27 [e] 

[O 2 ] [e] -^^31 +^34 



( 44 ) 



Nj* is produced by the following reaction: 

N^+hv - N2+e ( 45 ) 

The following reactions result in N 2 * loss: 



N 2 +O - NO*+N(^D) 


( 46 ) 


N 2 +O - 0*-^N2 


( 47 ) 
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( 48 ) 






Nz+e - N(^D) +N(^D) ( 49 ) 

Combining these reactions and solving for [N 2 *] yields the 
following solution: 



[Nz] 



^(^2) 

+K^-Jt,[0] [O2] +iCi5 [e] 



( 50 ) 



The following reactions lead to the production of O 2 *: 



Oj+Av - Oz+e 


( 51 ) 


n;^o^ - o;+jv2 


( 52 ) 


0**0^ - O 2 +O 


( 53 ) 



The following reactions lead to the loss of O 2 *: 



Oz-^NO - N0*-^02 


( 54 ) 


Oz+Ni^S) - NO* + 0 


( 55 ) 


Oz+e - 0+0 


( 56 ) 



Combining reaction (51) through (56) results in the following 
solution: 
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K 22 [NO] +JC23 [Ni^S) ] +it24 



(57) 



[o:i - 



The above equations are used in photochemical models 
such as the ones presented in Cleary (1986) . 

B. THE NIGHT AIRGLOW PHENOMENA 

The night airglow spectrum has been extensively studied. 
Far ultraviolet measurements of the night airglow spectra have 
been assembled from several sources including Apollo 16, 
numerous rocketborne experiments, and several space shuttle 
and Skylab based experiments. The nature of the spectrum is 
well established. For the purposes of this paper, the far 
ultraviolet night airglow spectrum can be divided into two 
regions; 1000 - I 6 O 0 A and 1600 - 2000 A. This serves two 
purposes; first it separates the spectrum into the respective 
regions of coverage for the two cameras, and second, it 
separates the region of primarily oxygen related emissions 
from the region of primarily NO related emissions. 

At night, the free ionospheric electrons and ions 
recombine resulting in a reduction of ion density. During 
this recombination process the atom or molecule can emit 
photons at various wavelengths. Figure (6) illustrates the 
night airglow spectra between 750 and 140oA. This region is 
characterized by the spectra of O, O*, He, and H. Of note are 
the atomic oxygen recombination emission lines at 1304A and 
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1356A. These are the only significant emissions in the 
filtered Camera One passband. 




Figure (7) shows the energy level diagram for 01. The 
transitions which are detectable by filtered Camera One cover 
the range from 10 to ‘7.5 eV. 

The nightglow spectra of the 1600-2000A region is 
characterized by the emissions of the NO-6 band. Figure (8) 
illustrates the nightglow spectrum from 1600 - 295oA. Note 
the presence of two NO 5-band emissions below 200oA. These 
are the only emissions detectable by filtered Camera Two, the 
upper limit of the camera prevents detection of NO y-band 
emissions and allows the study of 6-band emissions only. 
Figure (9) shows the energy level diagram for NO. The 
passband of filtered Camera Two equates to potentials between 
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Figure 7: Energy Level Diagram for 01. 

(Rees, Appendix 3) 



7.75 to 6.2eV. In Figure (9), the N and O atoms approach 
along the a^Hcurve to the state. The C^U state decays to 
the ground state emitting the 6 bands. The reaction is: 

N(*S) +0(^P) -‘NOia^U) -‘NOiC^ll) ^NO(X^U) +hv (6-bands) (58) 

The minima for the C and X states do not have the same 
internuclear separation. This leads to the emission of 
several bands in the v'=0 transition. (Rees, pp. 140-146) 
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Figure 8: Nightglow Spectra from 1600 - 

295oA. Obtained from the S3-4 
Satellite. (Rees ppl41) 



25 



Potential energy ( electron volts) 




Internuclear d»star>ce (A) 



Figure 9: Energy Level Diagram for NO, NO*, 
and NO'. (Rees, Appendix 3) 
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THE FAR UV CAMERAS EXPERIMENT 



The Naval Research Laboratory's Far UV Cameras Experiment 
was flown aboard STS-39 between April 28 and May 6 1991 as 
part of the Air Force Program 675 payload. The purpose of 
the Far UV cameras was to obtain images of the upper 
atmosphere, the shuttle environment, and celestial targets 
including stars and diffuse backgrounds. The wavelength band 
of the cameras was 105oA to 200oA. Two cameras were used to 
allow observations in two overlapping passbands; 1050 - 160oA, 
and 1230 - 200oA. Filters were used on some images to allow 
imaging in two distinct passbands; 1230-160oA, and 1650-2000A. 

A. THE FAR UV CAMERAS 

The electrographic Schmidt cameras developed at NRL are 
designed to provide for maximum sensitivity for detecting 
diffuse sources in the far UV range. This is accomplished 
through the use of (a) low f-number optics, (b) high 
efficiency opaque alkali halide photocathodes, and (c) 
microchannel intensification to permit detection of single 
photoelectron events. The objective was to develop from 
these far UV images a method by which an altitude profile of 
the upper atmosphere can be determined. Both cameras have 
aperture diameters of 75mm and a circular f ield-of-view of 20° 
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in diameter. The opaque alkali halide photocathodes have much 
higher quantum efficiencies in the vacuum ultraviolet than do 
common semitransparent photocathodes. As a result, the two 
cameras have extremely high overall detection efficiencies. 
The difference between the two cameras lies in their diffuse 
source detection capability and their operating wavelength 
region. 

1. Camera One 

Camera one covers the 1050 - 1600 Angstrom range. 
This camera uses a LiF Schmidt corrector plate and a KBr 
photocathode. The wavelength range can be limited to 1230- 
160oA by using a CaFj filter. This permits the removal of the 
strong Lyman-a emissions at 1216A from the image thereby 
preventing saturation of the image by this emission. Camera 
One does not have a microchannel plate (MCP) intensifier and 
therefore single photoelectron events are not detected by this 
camera. An MCP is not useful at this wavelength region 
because the intensity of Lyman-a nightglow is on the order of 
2000 Rayleighs and an MCP would saturate unfiltered images. 
The diffuse source sensitivity of camera one expressed as 
units of optical density (on the processed film) per 
KiloRayleigh-second is estimated at 0.04 d/kR-sec at a 
wavelength of 1216A. 
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Camera Two 



The second camera covers the 1230-200oA range. This 
camera uses a LiFj Schmidt corrector and a Csl photocathode. 
The range can be limited to 1650-200oA using a Si 02 filter. 
Photocathode detection efficiency versus wavelength is given 
in Figure (10) . An (MCP) is used in Camera Two. This is 
required in order to intensify the night airglow of 
wavelengths longer than 123oA where the typical intensity is 
less than 10 Rayleighs. The diffuse source sensitivity for 
camera two is estimated at 4.0 d/kR-sec at 1304A Camera Two 
can detect a IR emission of 1304A light at an average density 
of .1 using an exposure time of 25 seconds. The longest 
exposure time was 100 seconds, therefore an emission intensity 
of .25R would be detectable by camera two. The minimum 
intensity of the far UV diffuse galactic background is about 
2.5R when integrated over the camera two passband. The 
minimum intensity of the night 0* recombination airglow 
(1356A) is about IR. As a result of the MCP intensification, 
single photoelectron events are observable by camera two. 
(Carruthers, 1992) 
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Figure 10: Photocathode Detection Efficiency vs 

Wavelength. 



B. DATA REDUCTION 

Nearly 400 images were recorded by the two cameras during 
the STS-39 mission. After flight the images were scanned into 
digital format using a microdensitometer. This instrument 
digitizes values of optical density d-log{IJ I) where 1^, and 

I are the incident and transmitted intensities of the scanning 
beam. The images were stored at NPS on a Sun Microsystem 
computer in the Flexible Image Transport System (FITS) format. 
The tapes are manipulated using Interactive Data Language 
(IDL) . 
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DATA ANALYSIS 



A. FORMAT AND METHODS 

1. Shuttle Position Data Files 

The shuttle position data is contained on the 
Postflight Attitude Trajectory History (PATH) data tapes. The 
PATH data tapes used for this project are held at Naval 
Research Lab Washington DC, Space Science Division. The 
contents of the PATH data tapes are described in detail in 
NASA publication ICD-I-TOP-002 dated 15 January 1988. The 
required material from the tapes was obtained in November 
1992. All required PATH data for each exposure taken during 
the PF03 block of the FAR UV cameras experiment is contained 
in each record. A PATH data record is recorded for every 10 
seconds of flight except when transmission problems prevented 
reception. Although some breaks in data were experienced, the 
PATH tapes provide enough information to determine shuttle 
position and attitude accurately for each airglow image. A 
description of the data fields contained in each PATH record 
is contained in Appendix B. 

a. Positional Analysis 

The determination of where the center of each 
image is pointed is critical to this project. Specifically, 
the angle between the local vertical and the pointing 
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direction of the camera must be known. From this the 
observation zenith angle (OZA) can be determined for any point 
in the f ield-of-view. Therefore, the position and attitude of 
the orbiter must be determined. Two separate coordinate 
systems are used to define the position and attitude of the 
shuttle. The Aries-true-of-date system shown in Figure (11) , 
and the orbiter body axis coordinate system. The Aries-true- 
of-date system is based on aligning the +x axis with the right 
ascension of the true-of-date equinox of Aries. The xy-plane 
is aligned with the true-of-date celestial equator. The +z 
axis is upward through true north, and the y-axis completes a 
right hand system. Right ascension is measured 0-360 degrees 
from the +x axis positive in the direction of the +y axis. 
Declination is the angle between the xy-plane and the radial 
vector to the orbiter and measured positive toward the +z axis 
and negative toward the -z axis. The position of Aries used 
in this system has been updated since the 1958 standard was 
established (hence, true-of-date) . The SHPOS data field from 
the PATH data record gives the orbiter position in the Aries 
true-of-date coordinate system. Right ascension (RA) and 
declination (DEC) can be calculated from the xyz coordinates 
given in SHPOS. 

The second system is the orbiter body-axis system. 
This system is illustrated in Figure (12) . In this system the 
orbiter x-axis is aligned longitudinal to the body and is 
positive out the front of the orbiter. The z axis passes 
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through the center of gravity and is directed out the top and 
bottom of the orbiter normal to the x-axis. The z-axis is 
positive out the bottom of the orbiter. The shuttle y-axis is 
orthogonal to the xz-plane and completes a right hand system. 
The +y axis is directed outward toward the right wing. 




Figure 11: Aries True of Date Coordinate System. 

The SHXRA and SHXDEC data fields give the aimpoint of the 
shuttle +x axis in terms of right ascension and declination. 
The SHZRA and SHDEC data fields give the aimpoint of the 
shuttle -z axis in terms of right ascension and declination. 

Since the f ield-of-view is aligned with the 
orbiter -z axis, the orientation of the orbiter -z axis must 
be determined. The first step is to determine the orientation 
of the z-axis with respect to the orbiter radial (r) . Using 
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Figure 12 : Orbiter Body Axis 
Coordinate System. 



the unit vector defined by the RA and DEC of the orbiter, and 
the unit vector defined by the RA and DEC of the shuttle -z 
axis; the dot product of the two confirms the two vectors are 
orthogonal. This proves that the z-axis is orthogonal to the 
outward radial from the center of the earth. Using the 
Shuttle +x axis RA and DEC, the dot product of this vector and 
the radial confirm the two are parallel. Since both points 
pass through the shuttle position and are parallel the shuttle 
+x axis passes through the center of the earth. This verifies 
that the orbiter is nose down as expected. The final question 
is the position of the orbiter -z axis with respect to the 
Aries-true-of-date z-axis. To determine this, the RA and Dec 
of the orbiter, the RA and DEC of the orbiter -z axis, and the 
origin are used to define a plane containing the three points. 
The angle between this plane and the Aries-true-of-date +z 
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axis is the angle of the orbiter -z axis. Using a cross 
product of the two unit vectors to define a normal vector to 
the plane and the angle between the Aries-true-of-date +Z axis 
and the normal was calculated. This angle +/” 90 degrees 
gives the angle of the orbiter z-axis with respect to the 
Aries-true-of-date z-axis. Calculations confirmed the orbiter 
z-axis to be in the orbital plane which is at an inclination 
of 57 degrees. This confirms that the base of the shuttle was 
in the ram direction as expected. An illustration of orbiter 
attitude is given in Figure (13). Having confirmed the 
orientation of the f ield-of-view, the next step is to begin 
image processing. 
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(Carruthers, 1992) 
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2 . Image Processing 



Image processing is performed using Interactive Data 
Language (IDL) . The images were received from Naval Research 
Lab in the FITS format described above and were read using the 
IDL READFITS function. Four programs were developed. Three 
are required to recover data from an image. The fourth 
program is used to print output. The programs are contained 
in Appendix C. 

The first program, DVFILES.PRO, is activated by typing 
UVFILES, IMAGE at the IDL prompt. UVFILES.PRO is the procedure- 
name and the parameter IMAGE is the variable name of the 
image. This parameter may have any name as long as it is 
consistent for all four procedures. This procedure allows the 
user to access the image using the IDL PICKFILE widget. Once 
selected, the image is read using the READFITS function and is 
rebinned from 2048 x 2048 to 512 x 512 to facilitate 
processing. The image array is divided by 800 to convert from 
photodensitometer output to density units. A 512 x 512 window 
is opened and the image is presented to the user in black and 
white. The screen will prompt the user to select a color 
table. Color table (13) (Rainbow) is recommended. For 
exposures of 30 seconds, color table (1) will also give a good 
visual presentation of the image. 

Once the first program has run, the second program is 
called by typing UVBACKGROUND/ IMAGE. This program allows the 
user to remove background noise form the image. Noise can be 
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imparted to the image file by the PDS procedure described in 
the data reduction section, by copying the files over a data 
net, and background effects from the film that was used. The 
UVBACKGROUND procedure allows the user to define a box by 
selecting two points with the mouse cursor. This is done in 
one of the four corners outside of the circular image. Care 
must be taken not to include any portion of the circular image 
in the box and not to include any artifact or defect which 
will skew results. The mean density unit value of this region 
is calculated and this value is subtracted from the image 
array. The new image replaces the old image in the view 
window. 

The remaining two programs are called by typing either 
UVCURSOR, IMAGE or UVCURSORl ^ IMAGE . These programs produce 
plots of intensity verses observation zenith angle. The only 
difference between the two programs is that UVCURSORl produces 
postscript files of the output graphs for printing. Both 
procedures will produce on screen graphs. After calling the 
program the user uses the mouse cursor left button to select 
points on the image. The points will be connected by a line 
to aid the user. Up to 512 points may be selected. When 
completed, press the right button to exit the cursor loop. 
The program then performs calibration calculations. The 
program asks for the number of seconds of the exposure and for 
calibration factors in Rayleighs and Photons. Calibration 
factors are listed in Table (1) . The program then calculates 
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calibrated intensity values by multiplying the density units 
for each point selected by the cursor by the calibration 
factor and dividing by exposure time in seconds The UVCURSOR 
procedures produce plots of intensity versus observation 
zenith angle (OZA) . OZA is based on the aimpoint of the 
camera having an OZA of 90 degrees. The algorithm determines 
the number of degrees per pixel based on a 20 degree field-of- 
view and subtracts from a lower altitude limit of 100 degrees 
OZA. For future research, a correction factor to account for 
the camera aimpoint being not equal to 90 degrees can easily 
be added to the program. For the benefit of the user, Two 
plots are produced; the first gives intensity in Rayleighs, 
the second plot is in Photons/sec/steradian/cm^' The UVCURSOR 
procedure can be run 



Table I: FAR UV CAMERA CALIBRATION 
FACTORS. R - RAYLEIGHS, PH - PHOTONS. 



Camera 


Filter 


Factor 


1 


no 


3.596168E4 R 


1 


no 


2.858751E9 Ph 


1 


yes 


6.263172E4 R 


1 


yes 


4.978868E9 Ph 


2 


no 


7.074304E2 R 


2 


no 


6.623675E7 Ph 


2 


yes 


4.444141E3 R 


2 


yes 


3.532842E8 Ph 
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multiple times however the two plot windows must be quit 
first. After running the UVCURSORl procedure the graphics 
driver is set to postscript, this must be reset by the user to 
the default driver for his system before additional image work 
can be done. Therefore, the UVCURSORl procedure should be the 
last one used for a given image, and only once the user is 
ready for hardcopy plots. One additional plot is available; 
a contour plot. This plot is produced by the procedure 
UVCONTODR and is called by typing UVCONTOUR, IMAGE. This plot 
is in Rayleighs. As with UVCURSORl, this procedure should be 
run after UVBACKGROUND and the graphics driver must be reset 
after the procedure is run. 

B. RESULTS 

Six airglow images were obtained from NRL, all are filterd 
Camera Two images. The files are summarized in Table (II) . 

Table II: FAR UV FILENAMES, EXPOSURE LENGTHS, AND TIME OF 

IMAGING. 



Filename 


Exposure Length 


GMT 


FARUV462.IFITS 


3 seconds 


05.04.16.21.59 


FARUV143.IFITS 


10 seconds 


05.04.16.22.09 


FARUV143.IFITS 


30 seconds 


05.04.16.22.39 


FARUV62.IFITS 


3 seconds 


05.04.16.23.11 


FARUV63.IFITS 


10 seconds 


05.04.16.23.21 


FARUV55.IFITS 


30 seconds 


05.04.16.23. 51 
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A color display of theses images is produced using 
DVFILES, IMAGE. Figure (14) shows a contour plot of onr of 
these images (FARUV144 . IFITS) . This plot was generated using 
UVCONTOUR/ IMAGE. Next, the orientation of the orbiter's +x 
and -z axes were established. From this attitude information 
a method for determining intensity versus observation zenith 
angle was developed. A typical example of this is shown in 
Figure (15) . This figure shows the intensity profile 
corresponding to the rectangular area in Figure (14) . This 
profile and printout was produced using UVCURSOR/ IMAGE. 



41 



Contour Plot of For UV Airglow Intensity 
(Level X 100) 




Figure 14: Contour Plot of Intensity vs Bin. Rectangular 
Area shows the Vertical Column from which the Plot in 
Figure (15) was derived. 
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Intensity vs Zenith angle 




Intensity 

Rayieigns 



Figure 15: Intensity vs OZA for the Area shown in Figure 

(14) above. 
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V. 



CONCLUSIONS 



A. SUMMARY OF FINDINGS 

Shuttle position and orientation has been established for 
six images of the Far Ultraviolet Cameras Experiment. All 
images are filtered Camera Two images. The orientation of the 
orbiter's +x and -z axes were established and from this a plot 
of intensity versus observation zenith angle was produced. 

From these intensity plots, accurate measurements of 
thermospheric and ionospheric densities can be derived. From 
these densities an approximation of free electron density can 
be made. Since there are negligible numbers of negative ions; 
charge equilibrium requires that the number of positive ions 
equal the number of free electrons. Therefore, the density of 
positive ions equals the density of electrons. From this, a 
prediction of RF frequency effects can be made. 

B. TOPICS FOR FURTHER STUDY 

Since this is the first NPS thesis on the Far Ultraviolet 
Cameras Experiment, there are numerous topics for further 
research. 

1. Determination of thermospheric/ionospheric densities 
from the intensity plots and relating them to geographic 
regions. This would include the removal of background 
radiation from celestial sources. 
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2. Deriving an algorithm for processing the data from the 
Global Imaging Monitor of the Ionosphere satellite based upon 
the Far UV Camera images. Use of this algorithm to provide 
real time RF predictions to a geographic region. 

3. Spatial variations in charged particle densities in 
the ionosphere can be determined from spatial information in 
the images . 

4 . Any topics concerning the aurora using six additional 
images of the aurora obtained form NRL. This will require 
additional position analysis as all aurora images were 
gimbled. 

5. Additional unfiltered images are available for work 
concerning Lyman-a effects and 01. 
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APPENDIX A: PCE REACTION RATES 



Chemical Reaction Rates 

This section contains all the chemical reactions and reaction 
rate coefficients used for the density calculation programs 
for the neutral atmosphere model. Sources for the equations 
are given on the last page of the appendix. Equation numbers 
correspond to those in Cleary (1986) . 
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Chemical Reaction Rates (cont) 
Number Reaction Rate (cm^ s"’) 



1 


0 + hv 0*+ e 


I, = 9.1E-7 s*’ 


1 ^ 


1^2 -► 02* ^ 


I, = 1.6E-6 S’’ 


3 


Ng + hv - N 2 * + e 


I, = 1.2E-6 s’’ 


17a 


Nj* + 0 - NO* + N(^D) 


1.4E-10 (T/300)’°-*^ 


17b 


N 2 * + 0 - 0* + N 2 


l.OE-11 (T/300)’°*23 


18 


N 2 * + O 2 - O 2 * + N 2 


5.1E-11 (T/300)’°® 


19 


N 2 * + e -► N(^D) + N(^D) 


1.8E-7 (T/300)’°-^’ 


20 


0* + N 2 - NO* + N(^S) 


5.0E-13 


21 


0* + O 2 - 0 * + 0 


2.0E-11 (T/300)'°-'^ 


22 


O 2 * + NO - NO* + O 2 


4.4E-10 
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Chemical Reaction Rates (cont) 



Number 


Reaction 


Rate (cm^ s*’) 


1 


©2* + N('^S) - NO^ + 0 


1.8E-10 


24 


Oj* + e O + 0 


l.OE-5 (T)‘°-^ 


24 


NO* + e - 0 + N(^S) 


,25*4.2E“7 (T/300)*°*®5 


27 


NO* + e -► 0 + N(^D) 


.75*4.2E-7 (T/300)'°-®5 


29 


N(^S) + O 2 - NO + 0 


4.4E-12exp(-3220/T) 


29 


N(^D) + O 2 -» NO + 0 


6.0e-12 


30 


N(^D) + e - N(^S) + e 


6.0e-10 (T/300)®-® 


3* 


N(^D) + Q - N(^S) + C- 


5.0E-13 


33 


N(''S) + NO - N 2 + 0 


3.4E-11 


34 


N(^D) + NO - Nj + 0 


7.0E-11 
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Sources of Chemical Reaction Rates 



Equation number and source 

1. Torr et al (1979) 

2. Torr et al (1979) 

3. Torr et al (1979) 

17a. McFarland et al (1974) 

17b. McFarland et al (1974) 

18. Lindinger et al (1974) 

19. Mehr & Biondi (1969) 

20. Lindinger et al (1974) 

21. McFarland et al (1973) 

22. McFarland et al (1974) 

23. Goldan et al (1966) 

24. Torr et al (1976) 

26. Torr et al (1976) Kley et al (1977) 

27. Torr et al (1976) Kley et al (1977) 

28. Becker et al (1969) 

29. Lin & Kaufman (1971) 

30. Frederick & Rusch (1977) 

31. Rusch & Sharp (1981) 

33. Lee et al (1978) 

34. Black et al (1969) 
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APPENDIX B: PATH DATA FIELDS 

This appendix contains a description of the data fields 
used in the Post flight Attitude Trajectory History. A more 
complete description of the entire PATH product can be found 
in the PATH product documentation provided by NASA. 



Record number 


Record 


Orbit 


Orbit number 


GMT 


YYMMDD.HHMMSS 


ATTMAT 


3x3 Transformation 
matrix. 


SHPOS 


Shuttle position in the 
Aries true of date 
coordinate system. 


SHVEL 


Shuttle Velocity 
components 


SHXRA 


Shuttle X axis right 
ascension 


SHXDEC 


Shuttle X axis 
declination 


SHZRA 


Shuttle +Z axis right 
ascension 


SHZDEC 


Shuttle +Z axis 
declination 


SHLON 


Shuttle Longitude 


SHLAT 


Shuttle Latitude 


SHZSUN 


Angle between shuttle +Z 
axis and Sun 
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SHZEAR 


Angle between shuttle +Z 
axis and center of Earth 


L$SNVS 


0 - Sun not visible 

1 - Sun visible 


L$BBEC 


0 - Not in Electrical 
contamination band 

1 - In contamination band 


TANRAD 


Tangent radius 
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APPENDIX C: IMAGE PROCESSING PROGRAMS 



This appendix contains IDL programs used to read the Far 
UV images and create the plots of intensity versus observation 
zenith angle. 



pro uvf lies , image 
f ile=pickf ile (/read) 

result=readf its (file^head) 
image=rebin (result^ 512, 512) 
image=image/800 . 

window, 1, xsize=512, ysize=512 
tvscl , image 

read, (' select desired color table') 
read, clrtbl 
loadct , clrtbl 

return 

end 
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pro u vbackground, image 



cursor^ t , g, /device^ /down 
cursor, /device, /down 

j=min ( [ f , f 1 ] ) 
jl=max ( [f , f 1] ) 

k=min ( [g, gl ] ) 
kl=max ( [g, gl ] ) 

cornerl=image (j : jl,k:kl) 
areal=n_element s (cornerl) 
corner l=total (cornerl) 
cor nerl=cornerl /areal 



image=image -cornerl 
tvscl , image 

return 

end 
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} : • ! r , i ■ 

^ 1*‘M ivivii:;*' huttv'M lo \ int, j'less right button at laut point'} 



o u I ^ 1 , \ , y , /device, / dow' n 
couru or - 0 
a-fitarr (512) 
b^fltarr (512) 



while (!err ne 4) do begin 

a ( counter ) -y 
b (counter )= irrage ( X, y ) 

cur sor , xl , y 1 , /device, /down 
plots, [x,xl], [y,yl], /device 

counter-counter+1 
x*xl & y=^yl 

endwhile 



a=a ( 0 : counter ) 
b=b ( 0 : counter ) 

len=n_element s (a) -1 

degrees=20 . /512 . 
za=100-“ (a*degrees) 

read, ('Enter exposure tirie in seconds: '), seconds 

read, ('Enter Calibration Value for Rayleighs: '),calfactor 

read, ('Enter Calibration Value for Photons: '), photons 

intensity=b*calfact or /seconds 
intensityp=b* photons /seconds 

window, 0, xsize=512,ysize=512 
plot, intensity, za, psym^^l , $ 

xrange- [min (intensity), max ( intensity) ] , yrange* [max ( za ) , min ( za ) ] , $ 

title==^' Intensity vs Zenith angle', $ 

xtitle=' Intensity ! cRayleighs/sec/steradian/cm^2 ' , $ 

ytitle= ' Observation Zenith Angle in Degrees' 



window, 2,xsize=512,ysize=512 
plot , intensityp, za, psym=l , $ 

xrange* [min (intensityp) , max (intensityp) ] , yrange= [max(za) ,min(za) ] ,$ 

title=' Intensity vs Zenith angle', $ 

xtitle= ' Intensity ! cPhotons/ sec/steradian/cm'^2 ' , $ 

ytitle=' Observation Zenith Angle in Degrees' 



return 

end 
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V 



] ^ , { ' ' 1 ■ • ■ 1 ! t mou • button t (> s»* 1 f»ot point, press r i q h’ button at Iasi j .•*.•' t 



cursor, x,y, /device, /down 
cour. t ei = 0 
3=:itarr (512) 
b*fltarr (512) 



while t’.err ne A) do begin 

a (counter ) -y 
b ( count er )* image ( X, y ) 

cur sor, xl , y 1 , /device, /down 
plots ,[ X, xl 1 , [y, yl ] , /device 

counter-counter-*- 1 
x=xl & y*yl 

er.dwhile 



a=a(0:counter) 
b=t ( C ; counter ) 

len="n_element s (a) -1 

degrees=20 . /512 . 
za=lCC- (a ^degrees) 

read, ('Enter exposure time in seconds: '), seconds 

read, ('Enter Calibra^.ion Value for Rayleighs: '),calfactor 

read, ('Enter Calibration Value for Photons: '), photons 

intensity-b^calf actor /seconds 
intensityp=b* photons/ seconds 

window, 0, xsize=512, ysize=512 
plot, intensity, za,psym=l, $ 

xrange- [min(intensity) ,max(intensity) ],yrange-(max(za) ,min(za) ] ,S 
titles' Intensity vs Zenith angle' , S 

xt it le=' Intensity ! cRayleighs /sec/ steradian/cm‘'2' , $ 
ytitle=' Observation Zenith Angle in Degrees' 



window, 2, xsize=51 2, ysize=51 2 
plct, intensityp,za, psym-1 , $ 

X range- (min(intensityp),max(intensityp) ],yrange=[max(za),min(za) ],$ 

title-' Intensity vs Zenith angle' , $ 

xtitle-' Intensity ! cPhotons/ sec/ster adian/cm^2 ' , $ 

ytitle- ' Observation Zenith Angle in Degrees' 



55 



s»n y \ ' p3 ' 

d<.'v : vTo . : , I'^^nan'.e^ ' inlrnl . ps' , /port rail 



plot, :niensityp, 2 a, psym- 1 , $ 

Xiang e= (mm (intensityp) ,max (intensity p) ] , y range = [max ( za) ,min(3a) 
t It le*' Intensity vs Zenith angle', $ 
xtitle-' Intensity ! cPhotons/sec/steradian/ cm"2 ' , $ 
ytitle-' Observation Zenith Angle in Degrees' 

device, close 



devi re, : i lename=' inten2 .ps' , /portrait 
plot, intensity, za, psym=l , S 

X range = (min ( intensity ) , max ( intensity ) j , yr ange= (max ( za ) , min ( za ) j , 
title-' Intensity vs Zenith angle', $ 
xtitle-' Intensity !cRayleighs/sec/steradian/cm^2',$ 
yt i tie- ' Observation Zenith Angle in Degrees' 

device, close 
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p3 uvcont our , image 

r 3 d, 'enter exposure time in seconds' , seconds 
rad, 'enter Rayleigh calibration factor', cal 

iagel= image* cal /seconds 

iagel=smooth ( image 1 , 15) 



c t_plot, ' ps ' 

c’vice, f ilename=' uvcontour .ps 



' , /portrait 



C'ntour, imagel, leve 
t.tle=' Contour Plot 



ls=[50,100,150,200,250,300], c_annotation [ ' 
of Far UV Airglow Intensity ! c (Level x 100) 



. 5' 



civice , /close 



nturn 

nd 



1.5' , '2.0' , ' 
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